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ABSTRACT 



Context. Recent work suggested that the traditional picture of the corona above the quiet Sun being rooted in the 
magnetic concentrations of the chromospheric network alone is strongly questionable. 

Aims. Building on that previous study we explore the impact of magnetic configurations in the photosphere and the 
low corona on the magnetic connectivity from the network to the corona. Observational studies of this connectivity are 
often utihzing magnetic field extrapolations. However, it is open to which extent such extrapolations really represent 
the connectivity found on the Sun, as observations are not able to resolve all fine scale magnetic structures. The present 
numerical experiments aim at contributing to this question. 

Methods. We investigated random salt-and-pepper-type distributions of kilo-Gauss internetwork fiux elements carrying 
some lO^'^ to 10" Mx, which are hardly distinguishable by current observational techniques. These photospheric distri- 
butions are then extrapolated into the corona using different sets of boundary conditions at the bottom and the top. 
This allows us to investigate the fraction of network fiux which is connected to the corona, as weU as the locations of 
those coronal regions which are connected to the network patches. 

Results. We find that with current instrumentation one cannot reaUy determine from observations, which regions on 
the quiet Sun surface, i.e. in the network and internetwork, are connected to which parts of the corona through ex- 
trapolation techniques. Future spectro-polarimetric instruments, such as with Solar B or Gregor, will provide a higher 
sensitivity, and studies like the present one could help to estimate to which extent one can then pinpoint the connection 
from the chromosphere to the corona. 

Key words. Sun: magnetic fields — Sun; atmosphere — Sun: corona 



1. Introduction 

The appearance of the large-scale magnetic field on the 
quiet Sun photosphere is dominated by the so called net- 
work structure, which is found at the boundaries of the 
super-granular cells. This magnetic network with a typical 
scale of 20 Mm is made up by concentrations each carry- 
ing a magnetic flux in the range of 10^* to 10"*^^ Mx with 
field strengths of the order of kilo-Gauss l|Schriiver et all 
nWH) . The region in-between thes e network structures, 
the i nternetwork, is not field- free ijLivingston fc Harvevl 
Il975j) . but shows a small average flux density appearing 
to be of the order of a few to 50Mx/cm^. Depending on 
the spatial resolution, polarimetric sensitivity and diag- 
nostic technique (Hanle and Zeeman effect), authors fin d 
a weak volume filling component (Fa urobert et al.l|2001(l . 
or small concentrations with kilo-Gauss fiux tubes (e.g. 
Dominguez Cerdena et all l2003t Sanche z Almeida et al.l 



2004D . 



As there are fiux tubes with high, kilo-Gauss field 
strengths to be found in the network as well as the in- 
ternetwork, a better concept to distinguish between the 
network and internetwork would be the amount of fiux 
carried by each individual flux elements. With some 10^^ 
to 10^^ Mx each internetwork (flux tube) element car- 
ries about one to two orders of magnitude less mag- 
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netic flux than the network patches llLin fc Rimmeld 119991 : 
iSocas-Navarro fc Sanchez Almeidal l2002() . 

Even though each internetwork patch carries a very 
small amount of flux as compared to the stronger network 
patches, the weak internetwork patches can contribute a 
significant amount to the total (unsigned) fiux on the solar 
surface, because they are so abundant. In contrast to the 
old canopy concept (e.g. GiovaneUi, 1980), where the in- 
ternetwork patches have been neglected, the abundant low 
fiux but high field strength concentrations in the internet- 
work do interact with the strong network patches: a sizable 
fraction of the field lines originating from the internetwork 
patches do not close back to another internetwork patch 
very close by, but can connect either to the stronger net- 
work patches or to the corona. 

This was first pointed out by iSchriiver fc Titld l)2003|) 
who constructed a network patch surrounded by many 
smaller internetwork patches by placing magnetic charges 
at a plane, performed a potential field extrapolation and 
then studied the connections from their simulated photo- 
sphere, i.e. the distribution of point charges, to the corona. 
They found that, depending on the average absolute fiux 
density, only part of the network flux is actually connected 
to the corona. At an internetwork average absolute flux 
density of about 20Mx/cm^, which is compatible with so- 
lar observations, they found that only half of the network 
flux connects into the corona, the other half connects to the 
surrounding internetwork. Consequently, half of the mag- 
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netic flux in the corona originates from the internetwork, 
which is in strong contrast to the traditional canopy pic- 
ture. 

In this paper we will go one step further than 
ISchriiver fc Titlj l)2003(l and not only demonstrate that the 
coronal magnetic field to a large part originates in the in- 
ternetwork, but also investigate the role of the (spatial) 
distribution of the magnetic field in the low corona and the 
photosphere, i.e. at the upper and lower boundary, for the 
hnk between the photosphere and the corona. This investi- 
gation will show the limitations for the knowledge we might 
gain from magnetic field extrapolations on the connections 
from the network and internetwork to the corona. As it will 
turn out, this is basically through our imperfect ability to 
determine the photospheric magnetic field at small enough 
scales with sufficient polarimetric sensitivity. 

The organization of this paper is as follows. In order to 
also employ a boundary condition at the top of the com- 
putational domain we briefly describe a Fourier transform 
technique for the potential field extrapolation, the bound- 
ary conditions representing the photosphere and the low 
corona as well as how to derive the fraction of network 
flux reaching the corona (Sect.[2Jl. Based on these results 
for the extrapolations with different boundary conditions 
we first compare our results to previous work in Sect. 13. ll 
before we discuss our new findings on the impact of the 
low corona fSect. lX^ and the distribution of photospheric 
flux (Sect. lB^ on the magnetic connection from the photo- 
sphere to the corona. We will conclude the paper with a dis- 
cussion of the consequences of this study for further inves- 
tigations of the chromosphere-corona coupling in Sect.^ 



2. Magnetic field extrapolation and boundary 
conditions 

In order to determine the magnetic coupling between in- 
ternetwork fields and the lower corona we will calculate 
the fraction of the magnetic flux originating from the in- 
ternetwork photosphere reaching into the corona. In a pre- 
vious study ,Schrii ver 4^ Title (2003) computed "field lines 
by tracing the path of monopolar test particles through 
the summed potential field above a plane [resembling the 
photosphere] filled with magnetic point charges of mixed 
polarities that surround a small number of significantly 
stronger point charges" establishing a network like con- 
figuration. To enforce a vertical field at a coronal height, 
they added mirror charges at a height of 2=100 Mm. The 
go al of the present paper i s to re-visit the problem posed 
by ISchriiver fc Titla f2003) in order to further investigate 
the impact of the boundary conditions at the top, i.e. in 
the corona, and at the bottom, i.e. in the photosphere. 
Thus we will apply a magnetic field extrapolation allow- 
ing the allocation of reasonable boundary conditions at the 
bottom and the top of the co mputational domain. W ithin 
the framework of the study of lSchriiver fc Titlel ll2003 ) this 
would not be possible for the upper boundary. For the sake 
of simplicity we will be using a potential field extrapolation 
employing a Fourier transform to solve for the expansion of 
magnetic field into the half space a bove the photosphere , 
as has been described already by e.g. lAlissandraki^ l)1981|) . 



In a force free state, the Lorentz force has to vanish, i.e. 
VxB has to be parallel to the magnetic field. If one further 
assumes a current free state, then 

VxJ5==0. (1) 

implying that B can be written in terms of a gradient of 
a scalar field, i.e. S is a potential field. This has to be 
accompanied, of course, by 

V-B^O. (2) 

Starting from a given ^-component of the magnetic vec- 
tor in the z—0 plane and assuming a periodic nature of the 
problem in the horizontal directions, the field is calculated 
by utilizing a Fourier transform (FT) technique to solve Q 
and 121. 

2.1. Expansion of the magnetic field 

To solve the expansion of the magnetic field governed by 
^ and we follow Alissandrakis (1981), but in addition 
to his study we will also prescribe the (z-component of the) 
magnetic field at the upper boundary. 

The FTs, denoted by for the spatial derivatives 

of S = {BxT By, Bz) in the horizontal (x, y) and vertical 
(z) directions read 

J" [OA] = ik^T[B,] , T[dyB^] = ikyT[B,] , 
T[dA] = dzT[B.,] , 

and similar for By and B^- As the problem is assumed to 
be periodic in x and j/, the wavenumbers kx and fc^ are real. 
Applying these relations to the FTs of (^ and (j^l yields a 
second order homogenous differential equation 

dlBz = k^Bz . (3) 

Here the horizontal wavenumber is given through k — (fc^ -|- 
ky)^^"^ and for simplicity the hat symbol denotes the Fourier 

transform, i.e. 13^ = J^[Bz] - For the horizontal components 
we find 

B.^^dzBz , By^^dzBz. (4) 
The general solution of (|2Jl is 

Bzikx,ky,z) = ai e''^ + as e"''^ , (5) 

with ai=ai{kx,ky) and a2=a2{kx,ky). 

When considering a computat ional domain s tretch ing 
from zero to infinity (as done by lAlissandrakij. Il981|) in 
the vertical direction, naturally ai has to be zero! However, 
if we confine the considered domain to a certain height we 
can specify a lower and an upper boundary condition (BC). 
Based on the upper and lower BC we know the Fourier 
transforms of the vertical component of the magnetic field 
at the bottom and the top of our computational domain, 

= Bz{kx, ky, z—0) , 
Bzh = Bz{kx,ky,z^li) . 

This determines the coefficients in the general solution 

_ BzQ - Bz, „ _ Bz„ ~ Bz, e-^^ 
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All three components of the (FT of the) magnetic field can 
now be written in terms of the assigned boundary condi- 
tions and height z, 

^ ~ikx Bzf^ sinh[fcz] — B^q sinh[/c(z— ft,)] 
^ k smh[kh\ ' 

—iky B^^ sinh[fcz] — B^g sinh[fc(z— /i)] 
" smh[kh] ' 

^ Bz,^ smh[kz] — Bzg sinh[fc(z— ft,)] 
^ sinh[fcft] 

The magnetic field B is then determined by an inverse FT. 
2.2. Assigning Boundary conditions 

To achie ve a maximum degree o f comparability with the 
study of ISchriiver fc TitI3 l|200,'^ . we employed magnetic 
configurations at the lower boundary, z=0, which are de- 
rived by placing magnetic charges below the surface as 
was done in their work. From these we extracted the z- 
component of the magnetic vector at the bottom (z=0) of 
several computed field domains to create sets of magne- 
tograms for the use as lower BCs with the extrapolation 
model. 

2.2.1. Configuration of the lower boundary 

In order to establish network-like magnetic patches we as- 
sume a strong concentration of magnetic flux of 3-10^^ Mx 
in the middle of the horizontally periodic computational 
box, extending 14 Mm x 14 Mm horizontally. This is then 
basically the same setup as inlSchriiver & Title ( 2003), who 
considered 5x5 magnetic patches to achieve some period- 
icity. 

Surrounding the central patch about 200 smaller 
charges of equal size were randomly positioned to constitute 
the internetwork field. Each of these small patches carries 
the same absolute value of magnetic flux, but half of them 
has one, the other half the opposite magnetic polarity, i.e. 
the net flux from the background, viz. the internetwork, is 
zero. 

We generated the magnetic field distribution at the 
lower boundary by placing magnetic charges for the strong 
network patch as well as for the smaller patches just 
below the surface (internetwork z=— 124 km & network 
220 km) and mirror charges at -1-100 Mm. Then we 
could compute the magnetic field by summing up the mag- 
netic potential at z=0. In this way we ensured to have 
a lower boundary condition closely matching the study of 
[Schriivcr & Title (2003). This results in a horizontal exten- 
sion of the central network patch at the photospheric level 
(full width at half maximum, FWHM) of about 170 km. The 
surrounding weak internetwork patches have a FWHM of 
about 100 km. One can think of this lower BC for the mag- 
netic field as a magnetogram that should resemble a mag- 
netic network patch surrounded by a salt-and-pepper-like 
weak internetwork field. 

We constructed several of such synthetic magnetograms 
to simulate diverse flux distributions in the internetwork. 
They all have the same central network patch carrying 
3T0^* Mx, but they differ with respect to the spatial dis- 
tribution and the magnetic flux of the weak internetwork 



Table 1. Parameters of the lower boundary condition, 
i.e. for the strong network patch and the 200 surrounding 
smaller internetwork flux concentrations (see Sect. l2.2"T|) . 0; 
denotes the magnetic flux carried by each individual patch, 
A is the horizontal extent of the bottom boundary. 



network internetwork 



size of magnetic patches 
FWHM [km] 


169 




96 


flux of magnetic patches 
M [Mx] 


3-10i« 


10^5 . 


. .8-10" 


unsigned total flux 
El-^^l [Mx] 


3-10^® 


2-10" . 


. . 1.6-10^° 


unsigned mean flux density 
EI<^.|M[Mx/cm2] 


1.5 


0.1 


...80 


net flux 
E^' [Mx] 


3-10^** 








patches. We carried out experiments with the (absolute 
value of the) magnetic flux for each weak internetwork 
patch ranging from 10^^ Mx to 8T0^^ Mx. Please note that 
in each of the synthetic magnetograms all internetwork 
patches carry the same (absolute) flux. The internetwork 
thus has an unsigned mean flux density ranging from 0.1 
to 80Mx/cm^ (while the total flux is zero). The overall net 
flux density is always 1.5 Mx/cm^, since only the central 
magnetic concentration contributes to the net flux. (See 
TableP). 

Concerning the spatial distribution we constructed 22 
sets of randomly distributed internetwork patches. We then 
calculated a magnetogram for each of the 22 spatial dis- 
tributions with each of the 19 sets of internetwork fluxes. 
Together with the 4 conditions discussed in the next sub- 
section this adds up to a total of 1672 conflgurations. 

Through this we can study independently (1) the impact 
of the strength of the magnetic flux in the internetwork and 
(2) the role of the spatial distribution of magnetic patches 
in the internetwork, both at the bottom boundary, i.e. the 
photospheric level. 

2.2.2. Configuration of the upper boundary 

In contrast to the work of ISchriiver fc Titld ()2003(l we will 
also investigate the role of the upper boundary on the flux 
budget from the network patches into the corona. To ex- 
plore this we employed four different upper BCs for the 
magnetic field extrapolation outlined in Sect . 12. ll These are 
partly carried to the extreme (especially item 4) and are 
not necessarily realistic, but have been chosen to illustrate 
possible effects of the upper boundary. 

(1) In the first case we assume the field to be stretc hing 
to infinity (in the same way as e.g. I Alissandr aki^ ll 98 Ij) . i.e. 
ai=0 in (0. 

(2) The next case uses a Bz{x,y) configuration at 
2:=15M m extracted from the e xtrapolation method used 
also bv ISchriiver fc Titld ()2003D as an upper BC. By this 
we match the iJj-component at z=15Mm with their work 
for a reliable comparison. 

(3) In a third case we applied a 2D Gaussian curve on 
top of a constant field at a height of 2=15 Mm at the up- 
per boundary, i.e. a broad (FWHM 6 Mm wide) maximum 
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just above the central network patch. In order to have a 
reahstic Bz variation at the upper boundary, we examined 
an actually observed magnetogram (obtained with MDI) 
for a region that showed a magnetic configuration similar 
to what we have considered the quiet Sun in this work. 
Applying a standard potential field extrapolation, with an 
upper boundary as for case (1), we derived a variation of 
some w5.5% at 2;=15Mm. This we used for the amplitude 
of the Gaussian above the background. 

(4) For the last case we deployed a very unrealistic mag- 
netic configuration at z=15Mm to see what effect a really 
extreme upper boundary might have. For this we used the 
same Gaussian as in case (3), but now with —50% varia- 
tion, i.e. with the magnetic field strength being depressed 
by a factor of 2 at the upper boundary directly above the 
central network patch. 

2.3. Budget of magnetic flux 

Once three-dimensional grids of complete B{x,y,z) sets 
were computed through the extrapolation as described 
above we determined field lines by tracing paths of monopo- 
lar test particles through the box. 

2.3.1. Field line tracing 

iSchriiver &: Titld lj2003|) traced field fines starting m ran- 
dom directions from the central network patch in the z—0 
plane. Their method then follows magnetic field lines until 
they either connect back onto the source plane or reach a 
maximum height, i.e. an upper boundary. 

One of the main goals of our work is to find the fraction 
of magnetic flux from the strong network patch, that has 
made its way into the corona, i.e. to the upper boundary. 
Therefore we start our field line tracing at the top boundary 
(z=15 Mm). Once lines reach the z=0 plane, we count them 
either as connecting the network or the internetwork to the 
corona, depending whether or not they are coming down 
on the central network patch. 

We start tracing field lines from a regular grid at the 
upper boundary with a spacing of 100 km. By this reverse 
approach we achieve a high precision when answering the 
question where the source of coronal magnetic flux is lo- 
cated. 

For the purpose of plotting the filed line configurations 
as shown in Fig.^ and El we trace the field lines also from 
the bottom boundary, of course. Fig. ^ shows the field lines 
for a sample potential field extrapolation. At the bottom of 
the figure the distribution of the (vertical) magnetic field is 
displayed, while the upper part of the plot shows a mask, 
where the area at a height of 15 Mm connected to the cen- 
tral strong network magnetic field patch is plotted in white, 
while the area connected to the weak surrounding internet- 
work patches is shaded gray. 

2.3.2. Connecting network and corona 

The equivalent fiux carried by each field line counts either 
towards the network or the internetwork flux contingent. 
The amount of magnetic flux associated with a given field 
line is calculated through 

(pi = Bz{xi,yi,z^h) ■ Ai , 




Fig. 1. Composite plot of the magnetic field at the bottom 
(lower BC), field lines of the extrapolated magnetic field, 
and a mask of the area at 15 Mm height that is magnetically 
connected to the central strong network magnetic patch 
(white) . 



where Ai denotes the area of the resolution element at 
the top boundary to which the field line is connected to, 
i.e. the grid spacing of 100 km x 100 km. Respectively 
Bz{xi^yi,z—h) is the field strength at that resolution ele- 
ment. 

The sum of those fiuxes 0i associated with the field lines 
connecting the network patch to the upper boundary gives 
the flux $^gt^ in the corona originating from the magnetic 
network element. Comparing this to the total magnetic flux 
of the network patch <i>5vetwi i-^- 

/ = '^Srit^i^Nit^ , (6) 

yields the fraction of magnetic flux from the network which 
is connected to the corona. 

3. Results and Discussion 

We applied the extrapolation method and the calculation 
of the flux budget from the network patch into the corona 
for all of the 1672 sets of boundary conditions, viz. photo- 
spheric magnetog rams. As a subset these also include the 
approach used bv ISchriiver &: Titl3 l)2003j) . 

The fraction of magnetic flux that expands from the 
network patch to reach the corona was calculated for 
four distinct iJ^-configurations at the upper boundary (cf. 
Sect. 12.2.'^ combined with 22 different (random) distribu- 
tions of magnetic patches in the internetwork field (cf. 
Sect. l2.2.'T)) . Each of these 88 configurations was repeated 
for 19 different unsigned mean magnetic flux densities of 
the internetwork field ranging from 0.1 to 80 Mx/cm^ (cf. 
Sect. imil . 
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Fig. 2. Average fraction / of magnetic flux from the strong 
network patch into the corona as defined in © as a function 
of the unsigned mean internetwork flux density. The dot- 
da shed line shows the result for a point charge model similar 
to iSchriiver fc Titld l|2003|) . see Sect.O The thick solid 
line displays / for the field extrapolation using a Fourier 
transform technique. For all (!) of the four upper boundary 
conditions discussed in Sect. l2.2.'^ the results are practically 
the same. While the thick line represents the median value 
of / for the different spatial random distributions of inter- 
network patches, the thin lines show the scatter (enclosing 
2/3 of aU values for /). See Sect. lO and 1^ 



In the following we compare ou r results with those pre- 
sented bv ISchriiver &: Titlel l|2003|) in Sect.O investigate 
the impact of the upper boundary fSect. lT^ and study the 
role of the lower boundary fSect. lT5|) . 

3.1. Comparison to previous work 

Recalculating the models of ISchriiver fc Titld l|2003D in- 
cluding an evaluation of the flux budget enabled us to pos- 
itively verify their results. As they perform the extrapola- 
tion using a distribution of magnetic charges, their study 
might be called "point charge model". Based on such a 
model we calculated the average value of the fraction / 
defined in ^ of magnetic flux from the network patch be- 
ing connected to the corona for the various distributions 
of internetwork flux at the lower boundary fSect. imjl . 
The dot-dashed line in Fig.[21 shows the average fraction 
/ as a function of the unsigned mean internetwork flux 
density. We find the same roughly expon ential decrease of 
the fraction / as ISchriiver fc Titld l)2003 D. differing by 0.05 
to 0.08 as compared to their results. This deviation basi- 
cally stems from using either fleld line tracing starting from 
above (Sect. 12. 3.211 as d one in this work or the method of 
ISchriiver fc Titld l)2003(l starting from below. 

We can now also compare the point charge model with 
a potential fleld extrapolation utilizing a Fourier transfor- 
mation (FT) technique as outlined in Sect. l2.l1 with a up- 
per BC as given in case 2 of Sect. 12. 2."^ i.e. which has the 
same magnetic field distribution at the upper boundary as 
the point charge model. The result for the average frac- 
tion / for this experiment is plotted as a thick solid line 
in Fig.|5| The thin solid lines indicate the scatter of the / 
values derived for the different spatial distributions of the 




14.1 Mm 



Fig. 3. A selection of field lines starting from the inter- 
network are drawn for an unsigned mean internetwork flux 
density of 60 Mx/cm^. Here the computational domain is 
viewed from the side. The dashed line indicates the max- 
imum height of the highest-reaching closed field line from 
this selection. See Sect. l3.2l 



weak internetwork flux patches. The difference of the point 
charge model (dot-dashed) and the FT technique (solid) is 
rather small, which shows that the two methods give simi- 
lar results. Thus we regard it as justified that both methods 
are likewise suitable to tackle the problem at hand. In our 
work we proceed using the extrapolation model in order 
to be able to determine the possible impact of the upper 
boundary, which could not be handled by the point charge 
model. 

3.2. Impact of the lower corona on the network coupling 

We first turn to the impact of the coronal configuration on 
the connectivity from the chromosphere into the corona, 
i.e. the role of the upper boundary. Depending on the large- 
scale structures in the corona one might imagine different 
magnetic configurations at the low corona. For example, 
while one might expect a rather constant magnetic field at 
a height of some 15 Mm in a more or less unipolar coronal 
hole, one could expect a rather strong variation in a net- 
work patch located near an active region. To explore this 
we employ four different sets of upper BCs (see Sect. l2.2.2|) . 

For these cases we evaluated the fraction / of mag- 
netic flux from the network element reaching the corona 
in the same way as in the previous section. We find that for 
all (!) cases of the upper BC the resulting fractions / are 
practically the same. In other words, comparing the bud- 
get of magnetic flux connecting the network patch and the 
corona for substantially diverse iJ^-configurations at the 
upper boundary, we find no evidence of an impact of the 
magnetic field distribution in the low corona on the mag- 
netic connection from the chromosphere to the corona. This 
is true even for the quite "pathological" case 4 outlined in 
Sect.|2T3 with a very strong (unrealistic) depression of the 
magnetic field above the network patch by a factor of 2. 
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This result can be illustrated by a display of magnetic 
field lines when looking at the computational domain from 
the side (see Fig.|3)l- Above an altitude of about 5 Mm the 
magnetic field is basically open, i.e. all the field lines found 
above that magnetic transition are connected to the corona 
(the dashed line in Fig.O shows the height of the highest- 
reaching closed field line in that selection) . Above this mag- 
netic transition the small scale "salt and pepper" structure 
of the internetwork magnetic field is cleared and a more or 
less vertical i?-component dominates. The height of that 
magnetic transition is, of course, related to the mean dis- 
tance of the magnetic polarities in the photosphere, which 
is for the distributions of the present work as well as on the 
real Sun of the order of a couple of Mm. 

The good news of this little experiment is, that any 
physically reasonable assumption for the magnetic field at 
the upper boundary is sufficient, supposing it is assigned at 
an altitude above the magnetic transition, i.e. above some 
5 Mm. This is of special interest, as we are not really able 
to measure coronal ma gnetic fields as of yet, except for 
some s pecial cases fe.g. iLin et al., 1999t iLagg et all l2004t 
lWhitd .l200Fl. Thus even without a proper information on 
the magnetic field in the low corona we can study the mag- 
netic connectivity from the chromosphere into the corona 
by means of magnetic field extrapolations in a meaningful 
way on a granular and super-granular scale. 

3.3. Impact of the photospheric field distribution on the 
magnetic connectivity 

We will now turn to the role of the lower boundary, i.e. 
the photospheric magnetic field for the magnetic connec- 
tivity into the corona. As emphasized bv lSchriiver Sz Title! 
l)2003() a sizable fraction of the corona is connected not to 
the strong network patches, but to the weak internetwork 
fields. Thus when trying to relate coronal phenomena to 
chromospheric events, one cannot simply look for e.g. a 
correlation of the corona with a network patch below, but 
coronal events can also be triggered by processes in the 
internetwork. The basic question is here if we can use the 
magnetic field extrapolations to derive which coronal struc- 
tures are connected to which areas in the photosphere and 
chromosphere, no matter if in the network or internetwork. 

To tackle this question we will investigate the spatial 
distribution of the weak internetwork magnetic field patches 
in our simulated magnetograms, which are hardly resolved 
by curr ent instrumentation. For example MDI onboard 
SOHO llScherrer et all Il995|) has a sensitivity of some 20 
Gauss (if being optimistic) and a spatial resolution in the 
high-resolution modus of 1 arcsec («0.5 arcsec pixels) cor- 
responding to 725 km. Thus this instrument, which is often 
used to study the relations of magnetic field and coronal 
phenomena, can detect magnetic fiux concentrations down 
to some 10"'^^ Mx. Therefore most of the small internetwork 
flux patches ranging from 10^5 to S-lO^'^Mx as used in this 
study (cf. Sect . 12. 2. T]! are not detectable with MDI. And 
even for todays most sophisticated high-resolution spectro- 
polarimeters this is a difficult task, especially at the low flux 
end. For example POLIS at the German VTT in Tenerife 
lIBeck et al.Ll2nnl can go down to some 5 Gauss with a res- 
olution of 0.5 arcsec, if seeing conditions permit, i.e. down 
to fluxes of some 10^^ Mx. 

In conclusion, especially the small internetwork flux 
concentrations we inverstigatetd here are hard to resolve; 



with some 10^^ Mx they correspond to a 1000 G flux tube 
with 36 km diameter. As there is no real reason to be- 
live that such small structures do not exist, e.g. they are 
found in numerical simulation of magneto-convection (e.g. 
IVogler et al.Ll2005|) . it is an interesting task to investigate 
the influence of such small salt-and-pepper-like flux con- 
centrations to the magnetic connectivity into the corona. 

3.3.1. Connectivity and flux distribution 

To investigate the role of the actual spatial distribution 
of magnetic flux at the bottom, viz. the photosphere, it 
is instructive to study individual examples of the random 
distributions of the internetwork flux patches. 

Two examples are shown in Fig.^ The right column 
shows the lower boundary conditions for both cases, i.e. the 
magnetograms constructed from the random distributions 
of fiux patches of the internetwork and the network patch 
in the middle. The images display the vertical component 
of the magnetic field on a grey scale, with black and white 
representing opposite polarities. As outlined in Sect . 12. 3 .21 
we calculated the fraction / of flux from the network patch 
reaching the corona for various values of the flux of the 
internetwork patches (but respectively always for the same 
spatial distribution). 

The resulting fraction / as a function of unsigned mean 
internetwork fiux density is shown in the respective graphs 
in the left column of Fig.^l As can be easily seen the frac- 
tion of network flux reaching the corona is quite different 
in the two sample spatial distribution. While in one case 
(top row in Fig.^ basically all coronal flux is originating 
from the internetwork (!) for average unsigned internetwork 
fluxes above some 50Mx/cm^, in the other case always 
more than some 20% of the network flux reach the corona, 
even for very strong internetwork fields (bottom row). On 
average, i.e. when considering a large number of random 
spatial distributions of the internetwork patches, one ends 
up with a fraction / somewhere between the two extremes 
shown in Fig.^ and this average result was already dis- 
cussed in Sect. l3.ll and Fig.|5] Of course, this is simply due 
to the distribution of flux directly surrounding the central 
strong network patch. The middle column of Fig.^ shows 
the total internetwork flux surrounding the network patch 
in a square with a length of a side a, i.e. excluding the net- 
work patch. In one case (top row) the surroundings of the 
central (positive) network patch are dominated by nega- 
tive charges, in the other case (bottom row) the opposite is 
true. Thus depending on the dominating sign of magnetic 
fiux near the network patch one will get a very low or a 
very high fraction of the network fiux being connected to 
the corona. 

In a way, there is a chance for special configurations 
of the internetwork flux patches to "neutralize" the strong 
network patch, so that most of the coronal magnetic field 
actually originates from the internetwork. 

The intriguing part of this (in principle obvious) result 
is as follows. If the internetwork field patches are hardly 
resolvable by observations, but still can significantly alter 
the amount of fiux from the network reaching the corona, 
we are confronted with a clear limitation for current inves- 
tigations of the interaction between the chromosphere and 
the corona. 
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Fig. 4. Individual examples of spatial distributions of weak internetwork patches. Each row shows the results for one 
spatial sample distribution. Right colum: Lower boundary condition for the vertical component of the magnetic field 
shown as a magnetogram with black and white representing opposite polarities. The strong network patch is visible right 
in the middle of the magnetogram. Middle column: Amount of internetwork flux contained in a square with length of a 
side a centered around the central network patch for lengths a ranging from outside the central network patch to the full 
coverage of the horizontal extent of the computational box. Internetwork unsigned mean flux density: 40 Mx/cm^. The 
dashed line marks the magnetic flux of STO'^® Mx that emanates from the network element which in d is almost balanced 
by the internetwork in a short distance of roughly 2.5 Mm from the center. Left column: Fraction / of magnetic flux from 
the strong network patch into the corona as defined in (jSJ as a function of the unsigned mean internetwork flux density. 
Please note the significant difference in network flux connected to the corona for the examples shown here in the top and 
bottom row. See Sect. l3.3.Il 



3.3.2. Mapping the corona back to the photosphere 

So far we have examined the connectivity between the pho- 
tospheric field and the corona only quantitatively through 
the fraction / of network flux reaching the corona. We will 
now turn to mapping coronal regions back to the photo- 
sphere and vice versa. 

The masks in the left panels of Fig.|Sl show in white 
the areas in the corona at a height of 15 Mm which are 
connected to the central network patch in the photosphere, 
i.e. at the lower boundary. Here we used the same sample 
distributions of weak internetwork magnetic flux patches as 
in the previous subsection and Fig.^ 

This reveals that the corona is quite sensitive to the 
structure of the weak internetwork fields: the coronal vol- 
ume which is connected to the strong magnetic network 
patches changes significantly, if a different distribution of 
weak internetwork field patches is prescribed. Our compu- 
tations show that this is valid even down to very low flux 
concentrations (0.1-2-10"'^^ Mx) in the internetwork, which 
are hardly resolvable at the current state of observational 
capabilities. 

It is not only problematic to derive the coronal area(s) 
connected to the network patch, but also to identify which 
parts of the internetwork are connected to the corona seems 
hardly possible when using current instrumentation. This 



is illustrated also by Fig.El where the intersections with 
the top and bottom boundary of the fieldlines connecting 
the internetwork with the corona are indicated by circles. 
From the right panels it is evident, that the internetwork 
footpoints of the fieldlines connecting to the corona cluster 
at completely different regions, depending on the distribu- 
tion of the photospheric internetwork ffux. 

These results have serious implications for studies of 
the relation of the photospheric field to the coronal struc- 
tures, e.g. as been done in the context of solar wind ac- 
celeration an d its relat i on to the chromosph eric magnetic 
network fe.g. lXia et alll2n04lTri et alll2n05D . When using 
data from the MDI instrument, as mostly done in studies 
relating the corona to the photosphere in recent years, one is 
not able to resolve the weak internetwork fields as modeled 
in the present paper. This implies, that in those studies one 
cannot definitely pin down which areas in the low corona 
are in fact connected to which parts of the photosphere. 
Certainly, under the presence of weak fields, the connec- 
tion from the coronal features down to the photosphere 
will not be nice and funnel-type, but more spaghetti-type 
with a lot of connections from the coronal patches to the 
internetwork. Better future instrumentation with superior 
sensitivity for measuring photospheric fields will hopefully 
provide the necessary information to allow advanced stud- 
ies of the connectivity from the photosphere to the corona. 
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Fig. 5. Mapping from the photosphere to the corona. The 
left panels show in white the areas in the corona at an al- 
titude of 15 Mm which are connected to the strong cen- 
tral network patch for two distributions of internetwork 
magnetic patches (same as the examples shown in Fig.QJ. 
The black regions are connected to the internetwork. The 
right panels show maps of the vertical magnetic field at the 
bottom boundary, i.e. in the photosphere (re-plotted from 
Fig.EJ- Over-plotted are white circles which indicate the 
intersections of the magnetic field lines connecting the in- 
ternetwork patches with the corona at the upper and lower 
boundaries, i.e. the circles in the right panles show the 
footpoints of the fieldlines connecting weak internetwork 
patches with the corona. 



4. Conclusions 

In the present paper we investigated the role of the mag- 
netic field distribution in the photosphere and the low 
corona for the mapping from the chromospheric magnetic 
network into the corona utilzing magnetic field extrapola- 
tions. The classic picture was that the field fans out from 
strong network concentratio ns in the chrom osphere creat- 
ing an often called canopy ijGabriel 11976), which could 
be seen as a more or less stable found ation of the coro- 
nal m agnetism. Extending the work of ISchriiver & Titlj 
l)2003j) . who are challenging this scenario, we manipulate 
the boundary conditions of a potential field extrapolation 
at the top and bottom of a computational box stretching 
from the photosphere to the low corona at some 15 Mm 
height. This allows us to investigate the mapping from the 
photosphere into the corona in great detail. 

The good news are that the photosphere-corona con- 
nectivity is very stable against changes of the top bound- 
ary, i.e. the magnetic field in the lower boundary at heights 
of around 15 Mm. This is very reassuring when studying 
different regions on the Sun, such as coronal holes, active 
region patches or quiet Sun, where one might expect dif- 
ferent magnetic field structures in the low corona: when 
investigating the connection e.g. of coronal Doppler shifts 



and photospheric fields to study the driving of coronal dy- 
namics, a magnetic field extrapolation from an observed 
high-quality magnetogram seems to be sufficient to iden- 
tify the respective areas in the photosphere and the corona 
which are magnetically connected (cf. Sect. l3.2jl . 

However, the reliability of the magnetic connection de- 
rived from the field extrapolation hinges on the quality 
of the information about the weak small-scale magnetic 
field in the photosphere. Changes in the distribution of the 
weak internetwork magnetic flux patches, which are hardly 
observable with current instrumentation, can lead to dra- 
matic differences in the connectivity from the photosphere 
into the chromosphere, qualitatively as well as quantita- 
tively. Especially the mapping from the photosphere into 
the corona, i.e. which photospheric regions are magneti- 
cally connected to which coronal areas, is quite sensitive 
to the distribution of weak internetwork flux. With current 
instrumentation only beginning to resolve the structure of 
the internetwork magnetic field, one has to be careful when 
drawing conclusions on magnetic connections from network 
patches to the corona, as one cannot really tell about the 
role of the internetwork (cf. Sect. 13. 3(1 . 

It will be exciting to investigate the polarimetric studies 
with new high-resolution instruments from space as well 
as from the ground, e.g. the soon to be started Solar B 
mission or the 1.5 m solar telescope Gregor with their 
sensitive full Stokes polarimeters. Such studies will be of 
great importance to better understand the driving of the 
corona through the underlying photospheric convection. 
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